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ABSTRACT 

We present images from the Solar Blind Channel on HST that resolve hundreds of far ultraviolet 
(FUV) emitting stars in two kpc^ interarm regions of the grand-design spiral MlOl. The luminosity 
functions of these stars are compared with predicted distributions from simple star formation histories, 
and are best reproduced when the star formation rate has declined recently (past 10 — 50 Myr). This 
pattern is consistent with stars forming within spiral arms and then streaming into the interarm 
regions. We measure the diffuse FUV surface brightness after subtracting all of the detected stars, 
clusters and background galaxies. A residual flux is found for both regions which can be explained 
by a mix of stars below our detection limit and scattered FUV light. The amount of scattered light 
required is much larger for the region immediately adjacent to a spiral arm, a bright source of FUV 
photons. 


1. INTRODUCTION 

Numerous massive O and B-type stars, typically in 
stellar clusters and associations, form in the spiral arms 
of spiral galaxies. These stars emit strongly in the UV 
and thus strikingly delineate these spiral arms or armlets. 
However, UV emissi on is also present in t he inter arm re¬ 
gions of spirals fe.g.. lMarcum et ^1200111 . which is com¬ 
monly called “diffuse” emission because of its smooth 
appearance in Ultraviolet Imaging Telescope (UIT) or 
Galaxy Evolution Explorer (GALEX) images. Three dis¬ 
tinct phenomena likely contribute to this interarm UV 
flux. First, stars formed within the arms should disperse 
into the interarm regions due to their streaming motion 
with respect to the arms. Because O stars only live for 
^ 10 Myr, these dispersed stars should be predominantly 
B-type stars. Second, low-level star-formation may occur 
in situ within the interarm regions, directly contributing 
young UV-bright stars (the initial mass function [IMF] 
may or may not be the same as within the spiral arms). 
Third, dust within the interarm regions should scatter 
UV light from the strongly star-forming arms into the 
observed line-of-sight. 

UIT images originally showed that the FUV-NUV 
color of in terarm regions is redde r than that of the spi - 
ral arms (ILandsman et al.l Il992t ICornett et al.l 1199411 . 
The simplest explanation for this phenomenon is that 
the luminosity-weighted age of the stellar populations is 
older in the interarm regions than within the arms, as 
expected if spiral arms are the locations of increased star 
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formation. 

However, scattered light is also expected to contribute 
to the interarm UV flux. Monte Garlo-based radiative 
transfer models show that scattered flux in the FUV can 
range from approximately 0.1 to 1 times the measured 
stellar flux, depending on the type of dust, the dust ge¬ 
ometry, the dust dumpiness, and the optical depth to the 
midplane of stellar sources (iWitt fc GordonI [200011 . Ad¬ 
ditionally, a higher fraction of scattered light is expected 
when the viewing direction is close to perpendicular to 
the stellar disk, as is typically the case for measuring 
inter arm light. 

These model predictions fit with observations of the 
SMG and LMG which show that detected stars only 
make up 55% and 75% of the total FUV ligh t, respec¬ 
tively (iGornett et ahllTo^ iParker et al.l 119981 1. Detec¬ 
tion of linearly polarized light from the Wide-Field Imag¬ 
ing Survey Polarimeter (WISP) demonstrates that a sig¬ 
nificant portion of this unaccounted-for light is scattered 
(although exactly how much depends on dust-source 
relative distribu tion and the scattering phase function; 
IGo1e et alJlT^ . 


For UV photons to be scattered effectively, the opti¬ 
cal depth must be relatively low. As the optical depth 
increases, multiple dust interactions lead to a higher frac¬ 
tion of light absorbed. A low UV optical depth requires 
a low dust surface density, most likely fo und in the inter¬ 
arm o r outer parts of disk galaxies fe.g.. iHolwerda et al.l 
l2005f l or i n their halos, extending away fr o m the disk 
midp l ane ( Hodges- Kluck fc Bregmi^ I20l4 iSeon et HI 
12014 . iThilker et akf i 2005 ) found that the proportion of 
diffuse UV emission increased at large radii in the floc- 
culent spiral M33, as expected for scattered light. Ad¬ 
ditionally, diffuse UV emission is seen in the halos of 
the starb ursts NGC 253 and M82, far from the ir stel¬ 
lar disks (|Marcum et al.l[200ll : iHoones et al.ll200^ . This 
emission is most likely dust-scattered UV from their star- 
burst cores. 

Comparing GALEX images wi th data from t he In - 
frared Space Observatory (ISO), iPooescu et al.l (|2005f ) 
found that the far-infrared (FIR) to UV ratio is higher in 













































Figure 1. The two SBC pointings (small squares) over the 
GALEX FUV image of MlOl. 

the interarm regions than within spiral arms in MlOl and 
argue that this is best explained by dust-scattered light in 
the FUV. However, other explanations for a boosted in¬ 
terarm FIR to UV ratio are possible, including a shift to a 
slightly older stellar population. Here, we directly inves¬ 
tigate the scattering hypothesis to explain the interarm 
FUV light in MlOl. Our method relies on the excellent 
resolution of the Hubble Space Telescope’s (HST) Solar 
Blind Channel (SBC) detector, which resolves individual 
UV-bright stars at the 6.7 Mpc distance of MlOl. 

In Section 2, we describe the SBC data and our pho¬ 
tometric measurements of point-like sources from both 
the SBC data and optical Wide Field Camera (WFC) 
data. Section 3 presents our main analysis of the FUV- 
emitting stellar populations, including the recent star for¬ 
mation histories (SFHs) of our observed regions, color- 
magnitude diagrams and identification of stellar clusters. 
In Section 4, we evaluate the sources of the significant 
amount of remaining FUV emission in the two fields. 
Section 5 presents our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

2 . 1 . FUV Observations 

Using t he SBC de t ector of the Advanced Camera for 
Surveys (|Ford et al.l 119981) on the HST, we obtained 
deep, far-UV images with a 34.6" x 30.8" field-of-view 
of two interarm regions within MlOl. These two point¬ 
ings are shown in Fig. [1] The images were taken with 
the F150LP long-pass filter, which has an effective wave¬ 
length of 1614A and a FWHM of 177A. Pointing 1 was 
observed for a total of 10,200 s in a single visit of 16 ex¬ 
posures. Pointing 2 was observed for a total of 10,950 s in 
two visits of 9 and 8 exposures each. Small dithers were 
applied between individual exposures in order to reduce 
the impact of hot and warm pixels. 

While there is no read noise in the SBC detector, it 
does have a dark current that must be subtracted be¬ 
fore measuring the diffuse light. The detailed procedure 
we followed to remove the dark current from our data 
is described in Appendix A and summarized here. We 
find that the dark current can be modeled by two com¬ 
ponents, one that is always present (“low temperature”) 
and a second that appears only at higher detector tem¬ 
peratures (significant for approximately half of our data). 


Both components contain structure. The high tempera¬ 
ture stru cture peaks towards the center of the detector 
(see also iTeolitz et al.l (I2006D 1. while the low tempera¬ 
ture structure has a more modest gradient, increasing 
from the left to right hand side of the detector. We first 
create maps of the high temperature dark current, us¬ 
ing the SBC observations themselves, by subtracting an 
average and spatially smoothed low-temperature image 
from a similar high temperature image, leaving a map of 
only the high temperature component. This high tem¬ 
perature component is then removed from all appropri¬ 
ate individual frames, leaving only the cold temperature 
component. All frames are then combined together. The 
low temperature dark current is determined by compar¬ 
ing resolution matched SBC images to those of the same 
location taken with GALEX in the far-UV. 

After we remove the dark current, we combine all of 
the individual exposures together to make a single image 
of each field using DrizzlePac’s ASTRODRIZZLE with 
0.025" pixels. A WCS solution used to align the two dif¬ 
ferent sets of observations for Pointing 2 was computed 
using DrizzlePac’s TWEAKREG (STSci PyRAE pack¬ 
age). The resulting far-UV images are shown in Fig. [2j 

2.2. FUV photometry 

We identify point and point-like sources in the com¬ 
bined far-UV images using the DAOFIND task in IRAF. 
We use a detection limit of 4.5 counts, since the back¬ 
ground level is dominated by Poisson noise rather than 
read noise (observed to be below 0.5 counts per pixel). 
At the level of 4.5 counts, we are almost certainly detect¬ 
ing real sources {P < 0.02%). We found that we needed 
to run DAOFIND with two different sets of parameters, 
in order to detect all obvious sources. In the first run, 
the FWHM of the convolution kernel is set to 1.5 times 
that of the PSF (as commonly recommended for point- 
source detection in the read-noise limited case) and in 
the second to 2.5 times that of the PSF. Combining the 
two DAOFIND runs for each pointing results in initial 
source catalogs that include 709 objects in Pointing 1 
and 1523 objects in Pointing 2. 

We perform aperture photometry with DAOFIND’s 
PHOT task using a 4 pixel radius aperture, and use 
a zeroi^int of 20.4423 mag from the ACS Zeropoints 
websit^, in order to convert to the VEGAMAG pho¬ 
tometric system. We apply an aperture correction of 
—0.651 mag, determined from photometry of bright, iso¬ 
lated stars in the SBC/F150LP observations of globu¬ 
lar cluster NGC 6681 (Proposal ID = 11378). The ob¬ 
served luminosity function of the point (stellar) sources 
are shown as the data points in Fig. [31 The distributions 
increase steadily down to ^ 25.5 mag; at fainter magni¬ 
tudes than this, the source lists become incomplete. 

2.3. Optical Observations and Photometry 

Optical images of MlOl, covering the entirety of Point¬ 
ings 1 and 2, were taken with the Wide Field Camera on 
the ACS as part of Proposal ID 9490 (PI: Kunz). We 
obtained drizzled images in the following filters: F435W 
(- B), F555W (- V), and F814W (- I) from the Hub¬ 
ble Legacy Archive. The F435W and F150LP images 

® http://www.stsci.edu/hst/acs/analysis/zeropoints/#sbc 
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Figure 2. The SBC F150LP images of Pointing 1 (left) and Pointing 2 (right). 
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Figure 3. Distributions of F150LP apparent magnitudes of de¬ 
tected stellar sources (black points). The dashed black line repre¬ 
sents the expected F150LP magnitude distribution for a constant 
SFR (values given in Table 1). The grey band represents the range 
of random samples pulled from the best-fit distribution of SFH for 
each region. The vertical dotted line shows our estimated com¬ 
pleteness limit at 25.5 mag. 
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were aligned by identifying 20 bright stars in common, 
and using the IRAF task CCMAP to determine the off¬ 
set between the optical and far-UV images. We perform 


photometry on the optical images using a 3 pixel radius 
aperture for all sources brighter than mpi^ow = 25. Ze- 
ropoints of 25.791 (F435W), 25.738 (F555W), and 25.533 
(F814W) are used, as recommended by the ACS Ze- 
ropoint Calculator websit4I3- Aperture corrections of 
—0.28, —0.25, and —0.29 were applied to the F435W, 
F555W, and F814W filters respectively, based on the 
encircled energy dist ributions presented in the ACS In¬ 
strument Handbook (lUbeda fc et al.ir2012[) . We only in¬ 
clude optical photometry if the uncertainties are less than 
0.2 mag. Almost all sources brighter than mpibOLP = 23 
have good photometric measurements in all three opti¬ 
cal bands, while sources with tofisolp 25 only have 
good measurements 35%, 25% and 20% of the time in 
the F435W, F555W, and F814W images, respectively. 

3. ANALYSIS OF FUV EMITTING SOURCES 

In this section, we analyze the various sources of far- 
UV emission in our interarm MlOl fields: individual 
stars, stellar clusters, and diffusely distributed light. 

3.1. Stars 

All F150LP sources without an optical counterpart in 
F555W are considered to be stellar sources because any 
star cluster or galaxy contaminants would be detected 
in F555W. Specifically, the F555W images are sensitive 
to clusters down to 10^ Mq and background galaxies will 
have F150LP-F555W > 0 and thus be detected given the 
completeness limits. 

For F150LP sources with optical counterparts, we sep¬ 
arate stars from stellar clusters and background galaxies 
based on their size, since all but the most compact clus¬ 
ters are expected to be broader than the point spread 
function at the distance of MlOl (physical FWHM of 
3.2 pc). Source sizes were measured in the F555W im¬ 
age, because the signal-to-noise is generally higher in this 

http://www.stsci.edu/hst/acs/analysis/zeropoints/zpt.py 
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filter than in the FUV. Furthermore, mass segregation 
may cause the most massive stars (which emit strongly 
in the FUV) to sink to the centers of clusters, thereby 
leading compact clusters that have narrower FUV than 
optical sizes. In F555W, point sources should have a 
FWHM close to 2.0 pixels. Based on measurements of a 
hand-selected sample of stars, we find that point sources 
have measured FWHM < 2.35 pix. We adopt this simple 
FWHM approach because standard star-identification 
parameters (DAOPHOT’s “sharp” and “round”) do not 
work well in this extremely c ount-limited case. Further¬ 
more, iHolwerda et al.l (|2014f) demonstrate that source 
classification based on half-light radius (directly linked 
to FWHM for point sources) is more effective than other 
point-source identification techniques. Thus our source 
catalogs are dominated by individual stars rather than 
by stellar clusters. 

The shape of the stellar luminosity function in the FUV 
provides important constraints on the recent star forma¬ 
tion history. The observed luminosity functions of the 
stellar sources are shown as the data points in Fig.[3l with 
Poisson uncertainties. We have separated Pointing 1 into 
two distinct regions, since the north-west corner of this 
field contains the edge of a spiral arm with a markedly 
different luminosity function. The luminosity functions 
all increase down to the approximate completeness limit 
of 25.5 mag (shown as the dashed line in Fig. |3]). 

The shape of the stellar FUV luminosity function is 
driven by the recent star formation history. In order 
to interpret our observed distributions, we create pre¬ 
dicted luminosity functions that assume different star- 
formation histories. We populate these theoretical lumi¬ 
nosity functions star-by-star, first drawing the mass of 
the star from a Kroupa IMF, then assigning the age of 
the star according to the assumed star-formation history. 
We parametrize this recent star-formation history with 
a simple two-phase model with three parameters: the 
initial SFR, the final SFR, and a break time at which 
the SFR instantaneously changes from its initial to final 
value. This form of the SFH also provides an appropri¬ 
ate description for a scenario where young stars stream 
out of a star-forming arm into the inter-arm regions we 
have targeted. (In this case, the initial SFR maps to the 
within-arm SFR, the final SFR to a continuous level of 
SFR that may persist in the inter arm region itself and 
the break time to the time since the average FUV-bright 
star passed out of the spiral arm.) 

We use Geneva stellar evolutionary models to assign a 
F150LP magnitude to each star given its mass and age. 
The model isochrones provide log surface gravity {g), ef¬ 
fective temperature (Teff) and absolute V-band magni¬ 
tude (My) values for stars between 0.8 and 120 Mq at 
various agesE3- Synthesized sp ectra from Ku rucz stellar 
atmosphere models 1ATLAS9. iKurucj 1199311 then help 
us bridge the gap to a F150LP magnitude. As we are 
not concerned with absorption lines nor the ionizing con¬ 
tinuum, the Kurucz models should give nearly identical 
results to those of a non-LTE model such as TLUSTY 
for the B and early O stars (< 30 M(7i) of our sample (for 
comparisons and discussion, see iLanz fc Hubenvl (l2007f l 

We use the half-solar metallicity (z=0.008, 
ILeieune &: Schaerer! 1200111 mod els as approp riate based upon 
MlOl’s metallicity gradient (e.g.,|niOn|20TD. 


and iPrzvbilla et ahl (|2011IH . Given this equivalence, the 
Kurucz models were used for simplicity as they are in¬ 
cluded in the STSDAS synphot package. We select the 
appropriate stellar atmosphere spectrum for each star 
based upon its assigned log( 5 ) and Tgs and then com¬ 
pute the F150LP - V color by convolving this spectrum 
with Johnson V and FI50LP filters using synphot. With 
this color, we convert the My from the model isochrone 
into a F150LP magnitude. Binning these magnitudes, 
we obtain theoretical luminosity functions. 

We parametrize our recent star-formation histories 
with different break times (tbreak) of 6, 9, 12, 16, 20, 25, 
30, 40, 50, 60, 70, 85 Myr and ratios SFRi„itiai/SFRgnai 
of 0.01, 0.032, 0.1, 0.32, 3.2, 10, 32 and 100. Two addi¬ 
tional cases are considered: (1) A truncated recent SFH, 
with SFHflnai = 0 (which corresponds to no in-situ star 
formation in the inter-arm regions), and (2) a constant 
rate of recent star formation (which corresponds to a case 
where the SFR does not increase due to passage through 
an arm). We normalize SFRjnitiai such that the num¬ 
ber of predicted and observed sources match in the 25- 
25.5 magnitude bin. Table[T]lists the best-fit parameters 
for each location and each type of recent star formation 
history (broken, truncated, constant). The best fit was 
determined by minimizing the Poisson fit statistic: 

— 2 1nP = 2y^ TOi — rzi-I-Ui In—, (1) 

^ mi 

from lDolnhinl (1200211 . In this equation, we sum over each 
magnitude bin, where rii and mi represent the observed 
and predicted counts per bin, respectively. We fit using 
the magnitude bins up to 25.5. The grey bands in Fig. [3] 
represent the spread in fifty magnitude distributions ran¬ 
domly sampled from the best-fit parameters. 


Table 1 

Star formation history best fits 


Fit Type 

yi^'Hinitial 

(Mq yr-1) 

yFHfinal 
(Mq yr-l) 

^break 

(Myr) 

-21nP 

Pointing 1 - JMW 

Broken 

0.0073 

0.00073 

85 

4.8 

Truncated 

0.00076 

0 

6 

9.7 

Constant 

0.00073 

- 

- 

5.6 

Pointing 1 - SE 

Broken 

0.00157 

0.00016 

40 

0.5 

Truncated 

0.00084 

0 

12 

24.0 

Constant 

0.00069 

- 

- 

42.4 

Pointing 2 

Broken 

0.0037 

0.00037 

12 

3.5 

Truncated 

0.0036 

0 

9 

7.0 

Constant 

0.0032 

- 

- 

36.0 


Our star-formation history fits indicate that the two 
interarm regions (Pointing 1 - SE and Pointing 2) have 
very low amounts of ongoing star-formation consistent 
with SFRs being lower by a factor ft! 10 in the inter-arm 
regions when compared to those within an arm. Point¬ 
ing 2 is best fit by a star-formation history that breaks 
at 12 Myr, decreasing by tenfold at that time. Pointing 
1 - SE requires a an older break time of 40 Myr, but also 
a very low current SFR. The region corresponding to the 
edge of a spiral arm (Pointing 1 - NW) has the largest 
current SFR (despite covering the smallest area) and is 
well fit by a star-formation history that is either constant 
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or has been constant for about the past 85 Myr. Fit¬ 
ting with a Salpeter IMF instead of a Kroupa IMF does 
not change the best-fit SFHs, as expected because these 
two IMFs only differ below 0.5 Mq where stars do not 
emit appreciably in the FUV. We note that SFRs would 
be approximately 1.6 times higher if we had assumed a 
Salpeter IMF rather than a Kroupa IMF. 

Faint dust emi s sion from Spitzer and Herschel data 
(|Dale et al.l 120071 : iGordon et all l2008t iDale et ^120120 
prompts us to consider the effect of extinction on our UV 
luminosity functions. Dust surface densities in the two 
regions range from 0.5 to 2x10® Mq kpc“^ (Aniano et ah, 
in prep.), correspon ding to an Ay = 0.1- 0.4 by the em¬ 
pirical conversion of iKreckel et al.l (120131) . This conver¬ 
sion is valid for H ii regions, so based unon lCalzetti et all 
dHH, we expect a lower average extinction towards the 
stars {Ay = 0.06-0.26). We assume a typical value of 
Ay = 0.1 and conve rt this to ApisoLP = 0.26 using the 
iCardelli et al.l (I1989D curve with an Ry = 3.1. Applying 
this extinction to the theoretical models and refitting, 
the best fits for all three regions remain the same. 

The stellar color-magnitude diagrams (CMDs; Fig. 0]) 
qualitatively agree with the luminosity function derived 
SFHs, although they are limited in diagnostic power be¬ 
cause only the brightest stellar FUV sources are opti¬ 
cally detected (hence the far fewer number of stars plot¬ 
ted). The F150LP-F555W CMD (left panel) reveals that 
Pointing 1 contains more massive stars than Pointing 2, 
as expected based upon the star-forming spiral arm edge 
located in Pointing 1. Similarly, the recent ten-fold de¬ 
crease of star-formation in Pointing 2 can explain the 
greater number of evolved sources found in its CMD com¬ 
pared to Pointing 1. We note that the F150LP-F555W 
CMD (left panel) is much better at distinguishing masses 
for young, massive stars which are extremely degenerate 
in the F435W-F555W CMD (right panel). 

3.2. Clusters 

We have two cluster samples based upon two different 
goals. The first is designed to help account for all the 
FUV light detected and is thus very inclusive. It con¬ 
sists of all non-stellar F150LP sources other than the two 
sources in each field specifically identified as background 
galaxies. As a reminder, non-stellar F150LP sources were 
identified as those with optical counterparts with F555W 
FWHM above 2.35 pixels. Some of these sources may in¬ 
deed be background galaxies and some may be very low 
signal-to-noise point sources with inaccurately measured 
FWHM. The total amount of FUV light in these sources 
is 0.6 and 1.4% in Pointing 1 and 2, respectively. Thus 
this category of sources does not appear to contribute 
significantly to the interarm FUV emission in MlOl. 

The second sample is designed to investigate the prop¬ 
erties of the clusters themselves. For this sample, we 
identify stellar clusters from the F435W and F555W im¬ 
ages, since not all clusters are detected in the FUV. Clus¬ 
ter candidates are selected to have a FWHM larger than 
2.5 pixels in both filters, and to be brighter than 25.5 mag 
in the F555W band. This selection yields seven proba¬ 
ble clusters in Field 1 and ten in Field 2, eleven total of 
which are detected in the F150LP images. 

Fig. [5] shows color-color plots comparing the mea¬ 
sured colors of the clusters, denoted with diamonds 
or upper-limit arrows, with predictions from theoreti¬ 


cal evolutionary models. A very small correction for 
foreground extinction was applied to the observed col¬ 
ors (approximately one tenth the length of shown ex¬ 
tinction arrow). The solid (dotted) line represents the 
evolution of a Z=0.008 (Z =0.004) population as it ages 
(iBruzual &: Charlotll2003L hereafter BC03), with the log¬ 
arithmic age in years given above the line. The disk 
of M lOl has a steep metallicity gradient (e.g., iLi et al.l 
120131) . and these metallicities should approximate those 
at the location of our fields. The direction and length 
of the small arrow in the lower-left corner of each panel 
shows the expected change for an extinction of Ay = 0.3, 
an upper limit to the expected extinction. The top panel 
of Figure [5] plots optical colors (F435IU — F555IU vs. 
F555W — U8141U). While the data points are close to 
the model predictions, age and reddening are clearly de¬ 
generate in this color-color space. The lower panel, which 
includes the FUV measurement, helps to break the age- 
extinction degeneracy. Here, extinction will move clus¬ 
ter colors off the model tracks. However, in this panel, 
the cluster colors are reassuringly close to the predicted 
ones, indicating that there is very little extinction to¬ 
wards these objects. 

Estimating the ages of the clusters from Figure [SJ we 
find that they primarily formed with a range of ages be¬ 
tween 100 Myr to 1 Gyr ago. Only a single cluster has 
formed much more recently, within the past 10 Myr. The 
result that most of the clusters are of intermediate age is 
consistent with their location between spiral arms, since 
most clusters are expected to form within and not be¬ 
tween arms. These clusters have likely moved away from 
their birth sites due to their rotation with the MlOl disk. 

We estimate the mass of each cluster from the F555W 
luminosity and the predicted age-dependent mass-to- 
light ratio from the BC03 models. The clusters have 
a range between 10® up to a few times 10'^ Mq. The 
fairly low masses and intermediate ages of these clus¬ 
ters suggest that the outer MlOl disk is a fairly hos¬ 
pitable place for clusters to survive. Galactic shear will 
be minimal due to the flat rotation curve at these radii 
and tidal shocks from passing molecular clouds may be 
particularly infrequent. A full study of cluster masses 
and ages throughout MlOl would be required to more 
firmly establish the longevity of low-mass clusters in the 
outer disk. It is also possible that some clusters are even 
lower metallicity than we have considered here (if they 
are low-metallicity ancient clusters). In this case, their 
ages would be older and masses larger than we have cur¬ 
rently estimated. 

4. DIFFUSE LIGHT 
4.1. Estimates of non-detected stars 

The emission from non-detected stars is the most un¬ 
certain portion of this accounting. Main-sequence stars 
emit in the FUV up to about 1 Gyr and the large num¬ 
ber of these relatively lower-mass stars means they are a 
significant source of FUV light, despite their individual 
dimness. 

Our SFH fits only include information up to about 
100 Myr (Section 3.1). We estimate emission from non- 
detected stars in three different ways. The first way 
under-predicts the contribution from undetected stars 
by assuming that only stars younger than the directly- 


























6 




F150LP-F555W F435W-F555W 

Figure 4. Color-magnitude diagrams for the stellar FUV sources detected in the optical. Objects are color-coded by F150LP magnitude 
in both panels, with the scaling evident from the left panel. The dotted grey lines represent our completeness limits (mpisoLP = 25.5 and 
^F555W = '^F435W = 25.75). An extinction vector for Ay = 0.3 is shown as are error bars appropriate for each magnitude bin. Stellar 
tracks from z=0.008 Geneva models are plotted with the stellar mass of each noted in small font. Pointing 1 demonstrates a more extended 
upper main sequence than Pointing 2, while Pointing 2 exhibits a higher proportion of sources evolving off the main sequence. 


probed 100 Myr contribute to the FUV flux. The second 
case we consider our best estimate. It extrapolates the 
best fit SFRfinai (or the inter-arm SFR) back to 1 Gyr, 
under the assumption that the high SFRinitiai value is 
essentially a burst corresponding to passage through the 
spiral arm that does not last for the entire past 1 Gyr. 
(For Pointing 1 - NW which is best fit by a constant SFR, 
we assume a drop to one-tenth of the current SFR 100 
Myr ago for this case.) The third case extrapolates this 
higher SFRinitiai back to 1 Gyr and thus should be re¬ 
garded as an upper limit to the amount of non-detected 
FUV emission from stars. These three assumptions re¬ 
sult in the fractions graphically depicted for each region 
in Fig. 6 and tabulated for the best estimate (second 
method) in Table 2. 

As Fig. 6 shows, detected and non-detected stars can 
not account for all of the FUV emission in Pointing 1. 
The story is more ambiguous for Pointing 2. Pointing 2 
is a fairly isolated region, far from a spiral arm. A con¬ 
stant low level of star formation might be expected for 
such a region, however, there is a clear sign of a signifi¬ 
cant decline (by a factor of 10) in approximately the last 
10 Myr. Therefore, the best extrapolation of the SFR 
between 100 Myr and 1 Gyr is unclear, but the choices 
of SFRinitiai and SFR^nai likely bracket reality. In fact, 
using the higher SFRinitiai can be ruled out, because it 
would require more FUV emission than is observed. 

4.2. Diffuse emission as scattered spiral arm light 


Table 2 


Region 

PI - NW 

PI - SE 

P2 

Fa (F150LP) 

1.36E-11 

1.84E-11 

3.55E-11 

FUV diffuse SB 

1.76E-5 

1.57E-5 

0.00 

PAH SB 

0.29 

0.17 

0.09 

Fraction of FUV in: 




Background galaxies 

4.7E-5 

6.2E-3 

1.3E-3 

Glusters 

3.5E-5 

4.2E-3 

1.3E-2 

Detected Stars 

0.43 

0.059 

0.33 

Nondetected stars 

0.27 

0.32 

0.52 

Diffuse 

0.29 

0.61 

0.14 


For Av = 0.1 and SFRfinal extrapolated between 100 Myr and 
1 Gyr. Flux is in erg s“^ cm“^ A~^. FUV surface brightness is in 
counts s~^ pix”^; PAH surface brightness is in MJy sr~^. 

The radiative transfer models from iWitt fc Gordon! 
(|2000ll predict scattered to stellar flux ratios of 
Fscat/F* = 0.1 — 0.6 at an optical depth of ry = 0.5 for 
clumpy Milky Way dust at 1614 A (6.2 /rm see their 
Fig. 8). We note that the lower optical depth towards 
our regions (ry Ri 0.1) should give slightly lower scattered 
to stellar fractions, but the perpendicular viewing direc¬ 
tion should increase the same fraction due to increased 
probability of scattered light escape. Ratios of our mea¬ 
sured diffuse to stellar flux ratios for our preferred star 
formation history are 0.4, 1.6 and 0.16 for Pointing 1 - 
NW, Pointing 1 - SE and Pointing 2. While Pointing 
1 - NW and Pointing 2 are thus in the expected range. 
Pointing 1 - SE is higher. For this region, much of the 
scattered light probably does not originate within the 
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-0.5 0.0 0.5 1.0 1.5 

F555W-F814W 


Figure 5. Color-color diagrams for probable clusters in both 
fields. Lines show SSP evolutionary tracks at z=0.008 (solid) and 
0.004 (dotted) from BC03. The arrow in the lower-left corner shows 
the direction an extinction of j4v= 0.3 would move data points. The 
single digits mark logarithmic ages along the cluster evolutionary 
track. 


delineated region (instead from the spiral arm included 
in Pointing 1 - NW) explaining such a high fraction of 
scattered to local stellar light. 

We may also test for direct empirical evidence of scat¬ 
tered light in Pointing 1. As PI - SE is immediately 
adjacent to a UV-bright spiral arm, the amount of light 
scattered should depend both on the local radiation den¬ 
sity (and its direction) and on the density of scattering 
particles, here dust. 

We obtain dust su rface densitie s from Aniano et al. 
(in prep.) who apply [Praine fc Lil (I2007D models to the 
Spitzer and Herschel mid- to far-IR photometry of MlOl 
(along with several other galaxies in the KINGFISH sur¬ 
vey). A dust mass surface density of 1.2x10® Mq kpc“^ 
is typical for Pointing 1. Based on dat a for NGC 891 
(iBimchi &: Xilourii 120111) . NGG 4565 (|de Looze et al.l 
I2OI2D and s even other edge-on spira l galaxies with Her¬ 
schel data (|Verstappen et al.l I2013D , we adopt a scale 
height of 200 pc resulting in a dust volume density of 
3x10® Mq kpc“® in the disk midplane. Assuming uni- 
forml y distributed d ust and a dust cross-section per gram 
from IPrainel (|2003D , a FUV optical depth of r = 1 is 
reached within the midplane at only 210 pc, about a 
fifth of our SBG field-of-view. However, 200 pc above 
the midplane, r = 1 occurs after FUV photons travel 
about 560 pc. Thus, if dust is uniformly distributed, it 
is difficult for FUV photons to traverse much more than 
about 500 pc, even at the lower dust volume densities of 
the outer disk. 

Given these considerations, we look for a higher frac¬ 


Pointing 1 - NW 


SFH underestimate 


Detected Stai's 


Nondetected Stars Diffuse 


SFH best estimate 


SFH overestimate 


Pointing 1 - SE 


Pointing 2 


Figure 6. Showing the effect of different SFH assumptions on 
non-detected stellar flux. For each region, the top bar shows the 
minimal correction of accounting only for non-detected stars from 
the past 100 Myr (underestimate). The next two bars show the fit 
SFRfinal (best estimate) and SFRinitiai (overestimate) for the past 
100 Myr to 1 Gyr period. The small black region on the left-hand 
side depicts the fraction of FUV due to background galaxies and 
stellar clusters within the galaxy. 


tion of scattered emission in the region adjacent to the 
spiral arm in Pointing 1. To test this, we break the SE 
region of Pointing 1 into an intermediate region within 
210 pc of the spiral arm and a more distant region en¬ 
compassing the remainder of the field. Unfortunately we 
cannot refit for the SFH separately in these two regions, 
because so few stars are detected. Instead we assume 
the same correction factor for non-detected stellar flux 
in both. Performing this calculation, we find that the 
region adjacent to the spiral arm actually has less frac¬ 
tional diffuse emission and a lower diffuse FUV surface 
brightness. This is the opposite of the expected trend if 
scattered light dominates the diffuse FUV emission ob¬ 
served in our field. 

While this simple test fails to establish a clear signa¬ 
ture of scattered light, if the stellar population of the 
outer region is older than that of the inner region, the 
diffuse fraction in the above analysis would be under¬ 
estimated for the inner region and overestimated for the 
outer region, thus still leaving room for the expected pat¬ 
tern from scattering. 

4.3. Other sources of diffuse FUV 

Besides locally scattered FUV from MlOl, other poten¬ 
tial sources for the diffuse FUV emission include: fore¬ 
ground (Galactic) emission, background (extragalactic) 
emission, or emission from the warm ionized medium 
(WIM) within MlOl. 

Measurements of the full-sky diffuse FUV reveal that 



































































it correlates well with Hi column density and 100 
Jim, dust emission, implying a foregroun d origin (e.g., 
iParesce et alJll9^ iJakobsen et al1ll987ll . Based upon 
UV measurements from the Spectroscopy of Plasma 
Evolution from Astrophysic al Radiation (SPEAR) in¬ 
strument, ISeon et ^ (1201111 document the relation be¬ 
tween atomic hydrogen column density and /puv in their 
Fig. 21. Given the fore g round extinction predicted by the 
iSchlaflv fc Finkbeined (|2011ll model, we estimate that 
the foreground Hi column density is 4.3 x 10^® cm“^ 
in the direction of MlOl. At this low Hi column den¬ 
sity, the Ifuv curve gives a value of 350tioo photons s“^ 
cm“^ sr“^ A“^. Background galaxy emission is either 
in the form of identifiable background galaxies apparent 
in gaps of the GMC cloud cover (as the two galaxies in 
each field we have identified in the FU V) or is unlikely 
to hav e passed through the MlOl disk (iHolwerda et al.l 
l2007allbl l. Hence, we conservatively take 550 photons 
s“^ cm“^ sr“^ as an upper limit to the foreground 
and consider background FUV emission negligible. This 
value equates to 6.8 x 10“® counts s“^ pix“^ in our maps, 
approximately 500 times below the FUV surface bright¬ 
nesses that we measure. 

Within the Galaxy, the WIM contribution to the dif¬ 
fuse FUV from hydrogen t wo-photon emissi on is esti¬ 
mated to be about 4-9% (ISeon et al.l l20lHl . A pro- 
portionali t y wit h Ha emission was first put forth by 
iRevnoldi ((l992ll . here we use l 2 -f = 57.4 /na, where l 2 -f 
is in photons s“^ cm“^ sr“^ A“^, /hq, is in Rayleigh 
and the electron temperature assumed for the WIM is 
8000K. Using an Ha image of Ml Ol obtained with th e 
Kitt Peak National Observatory (jThilker et al.l 1200211 . 
we determine the average Ha surface brightness to be 
2.9 and 1.6 Rayleigh, for Pointing 1 and 2, respectively. 
Thus the expected contribution from two-photon WIM 
emission is 2.0 x 10“® and 1.1 x 10“® counts s“^ pix“^, 
respectively, also negligible. 

Recombination emission from Lya can pump molecu¬ 
lar hydrogen which then emits in the FUV band. This 
phenomenon is observed in planetary nebulae as well as 
young star forming regions. While the pumped-H 2 con¬ 
tribution to the diffuse F UV is around 10% in the Taurus- 
Perseus-Auriga complex (iLim et al.l 1201311 . we don’t ex¬ 
pect such a large contribution in our regions due to the 
lack of molecular hydrogen and the dimness of the ob¬ 
served Ha recombination emission. 

Having ruled out these other potential sources of dif¬ 
fuse FUV emission, we conclude that remainder of inter¬ 
arm FUV emission after stellar sources are removed is 
due to locally scattered light within MlOl. 

5. CONCLUSIONS 

Using deep and well-resolved FUV images from HST, 
we constrain the young stellar populations in two inter¬ 
arm pointings in the outskirts of MlOl’s disk. Hundreds 
of point sources are detected in both F150LP images, 
allowing us to create luminosity functions. These FUV 
luminosity functions have different shapes in the three re¬ 
gions studied. These shapes determine the best-fit recent 
star-formation history. In the two truly inter arm regions, 
we find evidence of truncated star-formation histories, 
indicating very little if any star-formation is presently 
occurring in these regions. 


Neither background galaxies nor stellar clusters signif¬ 
icantly contribute to the FUV light in interarm regions, 
amounting to a combined < 1.5%. However, several low 
mass (< 10^ Mq) clusters are detected in each pointing. 
Almost all have ages between 100 Myr - 1 Gyr, based on 
a comparison with stellar population tracks. The detec¬ 
tion of such clusters signals that the disruption of low- 
mass clusters is not extremely efficient in the outer disk 
of MlOl. 

In both interarm regions studied, significant amounts 
of FUV light remain after the detected stellar sources, 
clusters, and background galaxies are removed. The flux 
remaining is consistent with predictions of scattered light 
from radiative transfer models. In Pointing 1, which is 
immediately adjacent to a spiral arm, we find that scat¬ 
tered light from the arm must contribute significantly to 
the total FUV in that field (^ 60%). Meanwhile, Point¬ 
ing 2, which is further from any spiral arm or bright 
star-forming knots, has an upper bound of 16% scat¬ 
tered light. Thus the FUV emission in interarm regions 
in MlOl is partly due to individual FUV-bright stars and 
partly due to light scattered from the spiral arms. The 
latter contribution is only significant immediately adja¬ 
cent to spiral arms. 

The authors would like to thanks the referee for 
thoughtful comments which improved the manuscript. 
This paper is based on observations taken with the 
NASA/ESA Hubble Space Telescope obtained at the 
Space Telescope Science Institute, which is operated by 
AURA, Inc., under NASA contract NAS5-26555. These 
observations are associated with the program numbers 
9490 and 11147. 
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Figure Al. Dark current calibration images from 2008 organized by average detector temperature from left to right and top to bottom. 
The dark current shows a clear structure at higher temperatures, but the strength of the feature is not monotonic with average temperature. 
At lower temperatures variations are also present, most noticeably the globally higher value of the Avg. T = 8.48 observation. 


APPENDIX 

SBC DARK CURRENT REMOVAL 

In the ACS Data Handbook, the SBC dark current is listed as 0.8 —lx 10“® counts pixel“^. However, it also 
mentions that the dark current increases by about a factor of 5 over 2 hours (Section 4.4.1), as the detector temperature 
increases. SBC data obtained of the Hubble Ultra Deep Field also indicates this temperature-related in crease in dark 
curre nt, further noting that the dark current exhibits structure with a peak near the detector center (|Teplitz et al.l 

IM) 1 . 

We find that the dark current can be modeled by two compo nent s based on all the darks obtained with the SBC 
for 2008 (the year of our observations). These are shown in Fig. lAll having been binned to a 16 by 16 grid in order 
to increase the signal of the dark current. In this figure, the high-temperature dark images show a peaked structure 
that is absent in images taken when the detector temperature lower than about 9.5°C. Unfortunately, the detector 
temperature (as measured by the average of the two MDECODT keywords) does not appear to be the only parameter 
governing the dark current rate. Even at low temperatures, the dark current level can vary by a factor of two, although 
the structure of the dark-current in the low temperature darks appears similar, increasing from left to right. Based 
upon this inspection, we assume the dark current has a stable low-temperature structure (but arbitrary scaling) and 
an additional high-temperature component (with arbitrary structure and scaling). 

First we obtain high-temperature dark current maps by comparing our data taken at high temperatures (above 
9°C) to those at low temperatures. Fig. IA2I shows the detector temperatures of our sub-exposures as a function of 
time since the detector was turned on. The gaps show times when MIDI was unobservable during the orbit, so we 
average the two to five sub-exposures within these naturally-divided time blocks in order to increase the signal-to- 
noise. As the first time block of exposures is always quite cold (under 8°C), we assume these averaged frames have 
no high-temperature component dark current. For the higher-temperature averaged images (three such blocks for 
Pointing 1, two for Pointing 2 - Set 1 and only one for Pointing 1 - Set 2), we directly subtract the cold exposure 
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Figure A2. Detector temperature as a function of time since the SBC was turned on for our observations. Each exposure set has a similar 
profile with time and shows clear gaps where MlOl was presumably unobservable. Sub-exposure sets are grouped using these breaks to 
determine maps of the hot dark current. 


in order to obtain a map of the high-temperature dark current (always working on the coarse 16 x 16 grid, so slight 
positional uncertainties are not important). We then subtract these high-temperature dark current images from each 
sub-exposure within each corresponding high-temperature time block. All the frames for an observing day are then 
combined, producing an image with only the cold dark current component still present. 

We determine the scale of the cold dark current by comparing our SBC observations with dark-current free GALEX 
data (note the cold dark current structure is assumed to have the form of the averaged cold dark calibration frames). 
The SBC F150LP band is slightly redder than the GALEX FUV band (pivot wavelength of 1516A compared to 1612A). 
Because of this difference, we compute a color correction, based upon stars from the Bruzual & Chariot 1995 library, 
as provided by the STSdAs SYNPHOT package. We select stars that match the FUV-NUV colors we observe in the 
two fields, which range from spectral type A2 to B7. We use E(B-V)=0.0 and 0.5 and a Milky Way extinction curve, 
expecting these E(B-V) values to bracket the expected extinction in these outer regions of MlOl. The computed color 
correction is parametrized as: 


FxjSBC) 

Fx{FUV) 


1.66 - 0.594 


Fx{FUV) 

Fx{NUV) 


-h 0.129 


( FxjFUV) 

\FxiNUV) 


2 


(Al) 


Fx values are computed from provided PHOTFLAM values from iMorrissev et al.l (|2007f) and the ACS Zeropoints 
webpage. Typical color correction factors range from 1.03 to 1.10, with the maximum of 1.23. 

Applying this color correction to the GALEX FUV images leaves us with predicted dark-current free images for SBG 
F150LP at the GALEX resolution. Gonvolving the three SBC images with the GALEX FUV PSF using both the PSF 
and the IDL program of lAniano et al.l (|2011h and then resampling to the 1.5" pixels of the GALEX image, we may 
directly compare the SBC images to the GALEX prediction. To determine the dark c urre nt level, we fit a one-to-one 
line to the pixel flux density of the observed versus predicted image as seen in Fig. IA3I Such a one-to-one line is 
generally a good fit (some bright points off the relation are probably due to insufficient color correction for areas with 
bright OB stars), thus the offset provides the average cold dark current value. These offsets are 1.16 ± 0.04 x 10“^®, 
2.03 ± 0.02 X 10-1® and 2.20 ± 0.03 x lO'i® erg s"! cm'^ A"! in our three SBG exposure sets. These correspond 
to 0.74 X 10“®, 1.29 X 10“® and 1.40 x 10“® counts s”! in the detector pixels, comparable or slightly higher than 
the values given in the ACS handbook. We apply this scale to a normalized master cold dark computed from all of 
the low-temperature 2008 darks in order to remove the cold dark current from each of the sub-exposures. This dark 
subtraction accounts for approximately a quarter of the original flux in our images. 
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Figure A3. For the three different observing sets, the SBC measured flux density per l.h^'pixel compared to the predicted flux density 
from GALEX (using the FUV-NUV color correction to the FUV flux). Left to right: Pointing 1, Pointing 2 - Set 1 and Pointing 2 - Set 2. 





